AD-A066  935  INNSBRUCK  UNIV  (AUSTRIA)  INST  OF  EXPERIMENTAL  PHYSICS  F/G  A/2 

STUDY  OF  WATER  DRAINAGE  FROM  COLUMNS  OF  SNOW(U) 

JAN  1979  DENOTH,  A.j  SEIDENBUSCH,  W.  DA-ER0-76-G-038 

UNCLASSIFIED  PROJ.  4A161102AT24  TASK  A3  CRREL  79-1  N/L 


9-80 


66935 


I Study  of  water  drainage  from  columns  of  snow 


REPRODUCED  BY  ..... 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U S.  DEPARTMENI  OF  COMMERCE 
SPRINGFIELD.  VA.  22161 


Cover:  The  test  setup  for  drainage  of  snow  columns 
and  the  character  of  the  discharge  versus  time. 


CRREL  Report  79-1 


Study  of  water  drainage  from  columns  of  snow 


A.  Denoth,  W.  Seidenbusch,  M.  Blumthaler,  P.  Kirchlechner, 
W.  Ambach  and  S.C.  Colbeck 


January  1979 


UNITED  STATES  ARMY 
CORPS  OF  ENGINEERS 

COLD  REGIONS  RESEARCH  AND  ENGINEERING  LABORATORY 
HANOVER.  NEW  HAMPSHIRE,  U S A. 


Approved  for  public  release,  distribution  unlimited 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

I.  REFORT  NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NOJ 


CRREL  Report  79-1 


uL 


S.  RECIPIENT'S  CATALOG  NUMBER 


<8t 


TITLE  fvf 


S.  TYPE  of  report  a period  covered 


STUDY  OF  WATER  DRAINAGE  FROM  COLUMNS  OF  SNOW 


6.  PERFORMING  ORG.  REPORT  NUMBER 
S.  CONTRACT  OR  GRANT  HUMBERT*) 


fa 


AfDenoth,  W.jseidenbusch,  M.  blumthaler, 
P.jKirchlechnerJw.lAmbach  yd  S.G.  Gwlbteli 


(j^.  ^D^ERO-76-G-^38 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Institute  of  Experimental  Physics  -J  ( , 

University  of  Innsbruck  ^ 

Innsbruck,  Austria 


@di 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.S.  Army  Research  and  Development  Group,  Europe 
223-231  Old  Marylebone  Rd. 

London  NW15,  England 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 


DA  Proiect\4A161102AT24flask  A3, 
Work  Unit  001 


<31 


12.  BFPOHT  DATE 


Jam|hfy-W79 


z 


/ 


>1" NUMBER  8F  FWBPS 

18 


M.  MONITORING  AGENCY  NAME  A AODRESSfff  dll  lore nl  from  Controlling  Ottlco) 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
Hanover,  New  Hampshire  03755 


IS.  SECURITY  CLASS,  (ol  Me  roport) 


Unclassified 


ISa.  DECL  ASSI  FI  CATION/ DOWNGRADING 
SCHEOULE 


l«.  DISTRIBUTION  STATEMENT  (ol  Me  Roport) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  Ole  obmtroct  entered  In  Block  20,  II  dllloront  Irom  Report) 


IS.  SUPPLEMENTARY  NOTES 

Dr.  Colbeck’s  participation  was  funded  under  the  Directorate  of  Military  Programs,  Office,  Chief  of 
Engineers;  funded  by  DA  Project  4A161 102AT24,  Task  A3,  Work  Unit  001 


It.  KEY  VfOROS  (Continue  on  reeeree  elde  It  neeoeeery  end  Identity  by  block  number) 

Drainage 
Gravity  forces 
Snow 
Water 


. ABSTRACT  (Continue  on  reeeree  elde  II  neeoeeery  ai d Identify  by  block  ntmtber) 

Experiments  were  conducted  to  study  the  flow  of  water  through  columns  of  homogeneous,  repacked  snow.  The 
gravity  flow  theory  of  water  flow  through  snow  was  verified,  although  possibly  there  is  some  dependence  of  the 
relative  permeability  on  the  state  of  metamorphism  of  the  snow.  Also,  at  very  large  values  of  saturation  there  may 
be  some  additional  flow  in  saturated  channels. 

^ vw. 


^///22> 


oo  ,;sr„  M73 


EDITION  OF  t NOV  M IS  OBSOLETE 


TV 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Entered) 


PREFACE 


This  report  was  prepared  by  Dr.  A.  Denoth,  Dr.  W.  Seidenbusch,  Mr.  M.  Blum  thaler, 

Mr  P.  Kirchlechner  and  Univ.  Prof.  Dr.  W.  Ambach,  of  the  Institute  of  Experimental  Physics, 
University  of  Innsbruck,  Innsbruck,  Austria;  and  by  Dr.  S.C.  Colbeck,  Geophysicist,  of  the 
Snow  and  Ice  Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Eng.neer- 

' "^The'partkipation  of  Dr.  Denoth,  Dr.  Seidenbusch,  Mr  Blumthaler  Mr  Kirchlwhner 
and  Dr  Ambach  was  supported  by  a grant  of  U.S.  ARO  (Europe),  under  Pr0)ect  DAERO-76- 
G-038.  The  participation  of  Dr.  Colbeck  was  funded  by  DA  Project  4A161102AT24,  /?«e^c 
in  Snow,  Ice  and  Frozen  Ground-,  Task  A1 , Properties  of  Cold  Regions  Materials-,  Work  Unit 

001 , Properties  of  Snow  and  Ice.  . 

Dr.  G.D.  Ashton  and  Dr.  W.F.  St.  Lawrence  of  CRREL  reviewed  the  technical  content  of 


this  report. 


ii 


CONTENTS 

Page 


Preface ii 

Symbols  iv 

Introduction 1 

Test  procedure  1 

Theory  3 

Evaluation  of  the  experimental  results  5 

Example  evaluation 5 

Results  and  discussion 6 

Literature  cited ^ 

Appendix  A.  Measured  and  computed  discharge 9 


ILLUSTRATIONS 


Figure 

1.  Test  setup 2 

2.  Circuit  diagram  of  the  capacitative  level  pickup  and  the  threshold  switch  for  controlling 

the  magnetic  valve 2 

3.  The  theoretical  profile  of  discharge  rate  versus  time 4 

4.  The  theoretical  discharge  profile  of  total  discharge  versus  time 4 

5.  The  measured  water  discharge  points  and  calculated  profiles  of  discharge  versus  time 

for  various  values  of  n for  Experiment  no.  8 6 


TABLES 

Table 

1.  List  of  measured  and  computed  values 7 

i - 0.  % 


rj 


SYMBOLS 


a Abbreviation  for  AakS$n , m3/s  (rise  in  the  linear  part  of  the  V-t  diagram) 

A Cross  sectional  area  of  the  drainage  pipes,  m2 
d Grain  size,  mm 

g Acceleration  due  to  gravity,  m/sec2 

h Height  of  the  saturation  layer,  m 

k Intrinsic  permeability  of  snow,  m2 

L Length  of  the  drainage  pipe,  m 

n Exponent 

n Value  of  the  exponent  derived  by  interpolation 
P Fluid  pressure,  N/m2 

5 Licpid  saturation  (water  volume/pore  volume) 

S|  Irreducible  saturation 

5*  Effective  saturation,  S*  = (5-S,)/(1  -5,) 

t Time,  sec 

u Volume  flux  (water  volume/area • time),  m/sec 

Uq  Maximum  flux,  m/sec 

V Volume  of  the  escaping  water,  m3 

V„  Maximum  water  volume  which  can  escape,  m3 
Vs  Water  volume  in  the  saturation  layer,  m3 
Vc  Captured  water  volume,  m3 

W Liquid  water  content  (water  volume/total  volume) 

W s Irreducible  water  content 

z Vertical  coordinate,  measured  positively  downward  from  the  snow  surface,  m 
a Constant,  5.47  *106  m'1  sec’1 

ip  Porosity  (pore  volume/total  volume) 

<pt  <0(1-5,) 

p Density  of  wet  snow,  kg/m3 

pd  Dry  density  of  snow,  kg/m3 

pw  Density  of  water,  kg/m3 

p.  Density  of  ice,  kg/m3 

p Fluid  viscosity,  kg/m-sec 

Subscripts 
a Air 

c Capillary 

d Dry 

E End  of  experiment 

s Saturated 

v Delay 

w Water 

0 Initial  values 


iv 





STUDY  OF  WATER  DRAINAGE 
FROM  COLUMNS  OF  SNOW 

A.  Denoth,  W.  Seidenbusch,  M.  Blumthaler 
P.  Kirchlechner,  W.  Ambach  and  S.C.  Colbeck 


INTRODUCTION 

Previous  studies  of  the  flow  of  water  through 
snow  (e.g.,  Colbeck  1971 ) have  shown  that  gravity 
forces  dominate  the  flow,  at  least  in  well  metamor- 
phosed, ripe  snowpacks.  Colbeck  and  Davidson 
(1973)  found  that  water  flow  through  repacked, 
homogeneous  snow  could  be  characterized  by  a 
power-law  dependence  of  relative  permeability  on 
liquid  water  saturation;  in  particular,  the  relative  per- 
meability of  the  snow  to  the  liquid  increased  as  the 
liquid  saturation  rose  to  the  3.3  power.  For  simplicity, 
Colbeck  and  Davidson  (1973)  suggested  that  a third 
power  relationship  be  used  for  routing  meltwater 
through  snow. 

The  experiments  reported  here  were  carried  out 
with  snow  having  various  grain  sizes,  densities  and 
porosities  in  order  to  check  on  the  use  of  this  gravity 
flow  theory  under  different  conditions.  In  addition, 
an  attempt  was  made  to  determine  the  upper  and 
lower  limits  of  water  flow  rates  for  which  this  theory 
might  be  applicable.  At  very  low  flow  rates,  tension 
gradients  may  have  an  important  effect,  whereas  at 
very  high  flow  rates,  unstable  channelized  flow  may 
occur. 


TEST  PROCEDURE 

The  field  studies  were  carried  out  in  spring  on  the 
Seegrube  and  Hafelekar  areas  (1900  m above  sea  level 
or  2200  m above  sea  level),  near  Innsbruck,  Austria, 
and  in  the  summer  months  in  the  glacier  region  of  the 
Stubaier  Alps  (Schaufelferner/Daunferner,  3000  m 
above  sea  level),  during  the  period  April  1976  to 
September  1977.  The  percolation  of  meltwater  was 
studied  by  using  cylindrical  tubes  filled  with  sieved 
snow.  Initially  the  free  water  content  W0  and  density 


p0  of  the  snow  were  measured  and  the  grain  size  d0 
was  estimated  with  a lens  reticule.  After  being  filled 
with  the  sieved  snow,  the  tubes  were  set  upright  in 
deep  snow  to  guarantee  isothermal  conditions  during 
the  experiment. 

The  tops  of  the  tubes  were  closed  off  with  water- 
proof covers  so  that  it  would  be  possible  to  study  the 
process  of  drainage  of  the  free  water  already  present. 
The  tubes  were  covered  with  a snow  layer  0.5  to  0.75 
m thick  to  prevent  radiation-induced  melting  in  the 
snow  columns. 

Each  snow  column  was  supported  by  a wire  net  at 
the  bottom  and  a funnel  was  firmly  attached  to  the 
bottom.  After  several  days  of  drainage,  the  tubes 
were  weighed  to  determine  the  mean  density  of  the 
show  and  then  opened  to  allow  measurement  of  the 
water  content,  grain  size  and  density  of  the  snow  at 
various  heights  of  the  tube,  thus  permitting  the  deter- 
mination of  a complete  water  balance.  A sketch  of 
the  test  arrangement  is  given  in  Figure  1. 

The  water  content  was  measured  with  a device  de- 
veloped and  proven  in  the  Institute  of  Experimental 
Physics  of  the  University  of  Innsbruck.  This  device 
was  used  first  with  the  aid  of  a plate  condenser  to 
measure  the  dielectric  constant  of  the  snow;  the  free 
water  content  was  then  determined  from  a calibration 
curve  while  taking  into  account  the  snow  density  de- 
rived by  weighing  (Ambach  and  Denoth  1975). 

To  record  the  outflowing  water  voiume  precisely 
and  continuously  for  long  periods,  the  outflowing 
water  was  captured  in  a measuring  cylinder;  there  the 
water  level  was  measured  with  a pressure  probe  or  a 
capacitative  pickup.  The  measured  signal  was  used  to 
record  the  water  volume  (pen  record)  and  to  activate 
the  controls  of  a magnetic  valve.  In  this  way,  starting 
from  a preselected  height  to  which  the  snow  was  filled, 
it  was  possible  to  achieve  an  automatic  emptying  of 
the  measuring  cylinder  down  to  a preselected  lower 
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Figure  2b.  Circuit  diagram  of  the  threshold  switch 
for  controlling  the  magnetic  valve;  P2  is  the  lower 
switch  level  and  R is  the  relay  for  controlling  the 
magnetic  valve. 


level.  Figure  2 shows  the  circuit  diagram  of  the  con- 
trol electronics  with  the  threshold  circuitry  for  a 
capacitative  pickup.  A heater  was  provided  to  ensure 
constant  temperature  conditions  for  the  electronics. 

Because  the  device  used  for  liquid  content  meas- 
urements had  been  calibrated  in  the  laboratory  calor- 
metrically,  the  free  water  content  values  given  for  the 
individual  experiments  were  subject  to  an  absolute 
error  of  ± 0.5%  by  volume.  The  snow  density  p was 
determined  by  weighing  with  a relative  error  of  ± 1%. 


THEORY 

From  Darcy’s  law,  the  vertical  flow  of  air  and  water 
through  a porous  medium  can  be  described  by 
(Scheidegger  1957) 


•■■-as--) 


where  the  subscripts  "a"  and  “w”  stand  for  air  and 
water,  respectively.  In  these  equations,  ka  and  &w 
are  the  permeabilities  relative  to  the  fluids,  pa  and 
pw  are  the  fluid  viscosities,  Pa  and  Pw  are  the  fluid 
pressures,  z is  the  vertical  coordinate  measured  down- 
ward from  the  snow  surface,  and  g is  the  acceleration 
due  to  gravity.  For  the  vertical  percolation  of  melt- 
water through  snow,  these  equations  can  be  approxi- 
mated by  (Colbeck  1971) 

,3) 

where  Pc,  the  capillary  pressure,  is  normally  a func- 
tion of  liquid  saturation  and  porosity. 

Colbeck  (1971 ) wrote 

\ (4) 

pw\dS  dz  3 0 3z  f 

where  u is  now  taken  as  positive  in  the  downward 
direction  and  5 is  the  liquid  saturation.  The  first 
two  terms  are  generally  small  compared  with  the 
gravity  term  for  large  grained  materials  like  snow. 
Hence,  for  snow 


The  permeability  relative  to  the  water  phase  de- 
pends upon  the  intrinsic  permeability  of  the  snow/r 
and  the  liquid  saturation  S.  £w  approaches  zero  as 
liquid  saturation  approaches  its  irreducible  value  Sf 
and  kw  approaches  k as  the  pore  spaces  approach 
complete  liquid  saturation. 

Averjanov  (1950)  showed  that  for  porous  media 
6W  can  be  approximately  represented  as 


*w  = kS*n 


S*=(S-Si)/(]-Si)  (l 

and  S*  is  the  effective  saturation.  It  follows  that 
u = akS*n.  (< 


The  continuity  equation  is 

3"  + 0eM!  = o 

3 z 3 1 


0e  = 0(1-S,-)  (11) 

and  0 is  the  snow  porosity.  Combining  eqs  9 and  1 0 

n(ak)'ln  +0  >-’>/"  to  =0  (12) 

dz  dr 

which  describes  water  flow  through  snow  for  the 
case  where  the  snow  properties  do  not  change  signi- 
ficantly over  the  time  and  space  scales. 

The  solution  to  this  equation  can  be  given  in  the 
form  of  characteristic  slopes  where  the  speed  ( dz/dr ) 
at  specific  values  of  u is  (Colbeck  and  Davidson  1973) 

dS  =r(c‘k)'ln  u{nA)ln  (13) 

dt  u 0e 

where  d is  the  snow  grain  size.  For  decreasing  flux  as 
in  these  experiments,  it  follows  that  the  flux  at  any 
given  depth  is  given  by 


/0ez  \(n/*-1) 
\n(ak)'lnl 
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Figure  3.  The  theoretical  profile  of  discharge  rate  versus  time. 


Figure  4.  The  theoretical  discharge  profile  of  total  discharge  versus  time. 


where  t is  the  elapsed  time  since  that  value  of  flux  u 
left  the  surface  of  the  snow. 

At  the  start  of  each  experiment,  the  saturation 
throughout  the  snow  is  approximately  uniform  S0. 
The  initial  flux  of  water  is  w0;  thus  the  initial  dis- 
charge from  the  tube  of  snow  should  be  equal  to  u0 
for  a period  given  by 

f0  = — — 7 05) 

"MSS" 

where  Sfl  is  related  to  u by 


u0  = akStf’  (16) 

and  L is  the  length  of  drainage  pipe  (Fig.  1).  The 
theoretical  profile  of  u(t)  according  to  eq  14  is  shown 
in  Figure  3. 

The  volume  discharge  of  water  l/from  the  tube  is 
given  by 


i(“0t0+f  udt) 


where  A is  the  cross-sectional  area  of  the  tube.  It 
follows  that  the  total  discharge  is 


Vs  Zl|a/fSjf’r0+|-^_J'/(,'-,)(r7-1)  (18) 

(fl/(M_,1/{,-«>|  t>tQ 

This  theoretical  discharge  is  illustrated  in  Figure  4 
where,  in  the  limit, 

K.  = A Z.0(So-5j).  (19) 


EVALUATION  OF  THE  EXPERIMENTAL 
RESULTS 

One  object  of  this  data  analysis  is  to  determine 
the  best  value  of  the  exponent  n in  eq  7,  which  re- 
lates the  relative  permeability  to  the  liquid  and  the 
liquid  saturation.  The  sensitivity  of  this  exponent 
to  such  snow  parameters  as  grain  size,  density  and 
porosity  is  examined  as  is  the  general  nature  of  the 
usefulness  of  the  theory  of  flow  through  porous 
media  for  explaining  the  nature  of  water  flow 
through  snow. 

The  advantages  and  disadvantages  of  two  possible 
methods  of  evaluation  are  considered.  In  the  first 
method,  which  was  used  in  the  evaluation  of  the  ex- 
perimental data,  eq  18  is  rewritten 


V=  A UWQ-W-y(n-^{W^W^l(nA\  (20) 


AL(WQ-Wt 

* 0 

na 


a - AakSfi' . 


The  measured  volume  of  discharge  with  time  is  com- 
pared directly  with  the  theoretical  profiles  calculated 
from  eq  20,  where  A,  L,  W0  and  Wt  are  measured 
quantities  and  a represents  the  slope  of  the  measured 
discharge  curve  at  the  beginning  of  the  flow  ( t < t0). 
The  best  value  for  the  exponent  n,  the  only  unknown 
in  eq  20,  is  determined  by  interpolation  between  cal- 


culated profiles  for  two  selected  values  of  n,  bracket- 
ing all  of  the  measured  points.  The  advantages  of 
this  method  are  first  that  the  measured  volume  of  dis- 
charge is  used  directly  in  the  comparison  and  second 
that  the  permeability  k can  be  determined  from  eq  21 
with  high  accuracy. 

In  the  second  method,  used  by  Colbeck  and  David- 
son (1973),  eq  14  is  rewritten, 

\nu  = — + — Inf.  (22) 

"-1  n(ak)'ln  1-n 

The  use  of  eq  22  to  evaluate  the  experimental  data 
plotted  on  a log-log  scale  has  the  advantage  that  the 
important  snow  parameters  affect  only  the  intersect 
of  the  plot  whereas  the  exponent/;  can  be  deter- 
mined solely  from  the  slope  of  the  plot.  But  the  dis- 
advantage of  using  eq  22  is  that  the  flux  of  water  u 
must  be  determined  from  the  slope  of  the  discharge 
curve,  a procedure  which  may  lead  to  large  errors. 
However,  it  is  necessary  to  measure  W0  and  W-t  to  cal- 
culate n with  the  first  method  whereas  these  data  are 
not  used  in  the  evaluation  according  to  the  second 
method.  In  the  current  set  of  experiments,  the  first 
method  is  preferred  in  order  to  avoid  the  calculation 
of  the  rate  of  discharge. 


EXAMPLE  EVALUATION 

Of  the  10  experiments  in  this  study,  experiment 
no.  8 is  used  as  an  example  in  this  data  analysis  because 
it  is  typical  of  the  other  sets  of  experiments.  In  all 
of  the  experiments,  the  outflow  was  initially  delayed 
because  of  the  formation  of  a saturated  layer  at  the 
base  of  the  snow  column.  The  volume  of  the  satur- 
ated layer  is  given  by 


Vs  = atv 


where  a is  the  slope  of  the  discharge  curve  initially 
and  fv  is  the  delay  time  before  the  onset  of  discharge. 
Because  of  the  limited  capillary  effect  in  snow,  it  is 
assumed  that  the  formation  of  the  saturated  layer 
did  not  affect  the  flow  in  the  remainder  of  the 
column;  hence,  the  initial  flux  u0  prevailed  through- 
out the  column  at  first.  The  water  volume  Vs  of  the 
saturated  layer  was  added  to  the  measured  discharge. 

The  measured  profile  of  water  discharge  versus 
time  for  experiment  no.  8 is  shown  in  Figure  5.  Pro- 
files calculated  from  eq  20  for  n = 2 and  n = 3 are 
shown  to  bracket  all  of  the  measured  points.  By 
interpolation,  the  best  value  of  n is  found  to  be  2.38. 


Figure  5.  The  measured  water  discharge  points  and  calculated  profiles  of  discharge  versus 
time  for  various  values  of  n for  Experiment  no.  8. 


RESULTS  AND  DISCUSSION 

All  of  the  experimental  results  are  given  graphically 
in  Appendix  A in  the  form  of  measured  volume  of 
discharge  versus  time.  For  each  experiment,  cal- 
culated profiles  of  discharge  versus  time  are  shown 
for  different  values  of  the  exponent  n.  A compila- 
tion of  all  of  the  determined  quantities  of  n and  k by 
use  of  the  first  method  is  given  in  Table  I.  The  values 
of  n were  found  to  range  from  a low  of  2.16  to  a 
high  of  4.59. 

There  appears  to  be  a dependence  of  n on  the 
grain  sizes  or  the  sizes  of  clusters  of  grains  and  hence 
upon  the  age  and  history  of  wetness  of  the  snow.  For 
example,  experiments  no.  1 and  7,  where  large  clus- 
ters of  crystals  were  observed,  are  described  by  n 
values  of  4.59  and  4.31,  respectively,  whereas  for  the 
relatively  young  snow  of  experiment  no.  1 0,  the 
grain  size  was  about  0.5  mm  and  n was  2.16.  How- 
ever, no  clear  relationship  exists  between  any  snow 
parameter  and  n. 

For  experiments  no.  10  and  1 1 , the  initial  flow  of 
water  was  larger  than  expected  theoretically.  In 
these  experiments  the  comparatively  large  value  of 
initial  saturation  (~  1 5%)  could  have  caused  the  for- 
mation of  unstable  fingers  of  saturated  flow  in 
addition  to  the  unsaturated  flow  described  by  the  un- 
saturated flow  theory  given  here.  This  value  of  initial 
saturation  corresponds  closely  to  the  change  from  the 


pendular  to  the  funicular  mode  of  saturation  where 
some  of  the  assumptions  used  in  the  derivation  of  eq 
9 might  be  invalid.  It  is  very  important  to  note  that 
the  large  values  of  water  flux  associated  with  these 
large  values  of  saturation  are  greater  than  anything 
normally  encountered  in  nature.  The  question  of  un- 
stable flow  arising  from  high  saturations  may  be  im- 
portant in  experiments  where  large  flow  rates  are 
often  encountered  although  it  may  not  be  an  import- 
ant question  in  naturally-occurring  snowpacks. 

As  shown  by  the  experimental  data  and  computed 
curves  of  volume  versus  time  (App.  A),  the  gravity 
flow  theory  offers  a reasonable  explanation  of  the 
percolation  of  water  through  snow  for  the  wide  range 
of  flow  rates  and  snow  conditions  studied.  Also,  the 
values  of  snow  permeability  k calculated  from  eq  21 
are  compared  with  the  values  calculated  from  Shimizu’s 
(1970)  formula, 

/?  = O.O77</2exp(-7.8(1-0)pj/pJ  (24) 

where  the  mean  grain  size  and  the  porosity  are  shown 
in  Table  I.  These  two  methods  result  in  generally 
good  agreement,  suggesting  that  Shimizu’s  formula 
can  be  extended  to  a grain  size  of  about  3 mm. 

Further  investigations  of  these  phenomena  should 
concentrate  on  developing  an  understanding  of  the 
unstable  flow  of  saturated  fingers  and  the  mode  of 
water  flow  through  fresh  snow.  In  both  cases,  the 
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Table  1.  List  of  measured  and  computed  values. 


Experi- 

ment 

w0  “0 
(%)  (kglm3) 

<t>0  s0 
(st)  (Si) 

d 

(mm) 

A 

(m2) 

L 

(m) 

hs 

(m) 

(m3) 

(m3) 

1 

5.9 

500 

51.9 

11.4 



0.0625 

1.17 

7.7X10'* 

2.5  XI  O'4 

1.9X10'3 

2 

2.8 

500 

48.5 

5.77 

1-2 

0.0625 

1.12 

2X1 0'2 

5.5  XI 0'5 

5.976X10"4 

3 

8.7 

520 

52.8 

16.5 

1.5 

0.073 

1.4 

1.21X10'2 

4.67X1 04 

4382X1 0'3 

4 

7.0 

460 

573 

12.2 

0.5 

0.073 

1.3 

1.67  XI 0'2 

7X104 

238X1 0'3 

5 

7.0 

460 

573 

12.2 

0.5 

0.073 

1.4 

1.05  XI 0'2 

4.4X1  O'4 

3.05  XI 0'3 

6 

8.7 

560 

48.4 

17.97 

1-13 

0.073 

1.4 

3.06X1 0'2 

1. Of  XI  o'3 

3.78X1 0'3 

7 

8.0 

550 

48.7 

16.4 

1-3 

0.073 

1.4 

2.05  XI 0'2 

7.3X1  O'4 

3.36X1 0'3 

8 

7.8 

550 

483 

16.1 

1 

0.0625 

1.4 

1X10'2 

3X 10"4 

4.1 2X1 0'3 

9 

4.7 

540 

46.2 

10.2 

1 

0.0625 

1.34 

1.04X1 0‘2 

3X10"4 

1.57X1 0'3 

to 

7.7 

524 

51.3 

15.02 

0.5 

0.0625 

1.4 

2.5  XI 0'2 

8X 10"4 

2.3X1 0'3 

11 

7.3 

528 

50.4 

143 

1 

0.073 

1.4 

6X1 0'3 

2.2X1 04 

3.091X1 0'3 

Experi- 

ment 

Duration 

(s) 

vt*ve 

(m3) 

PE 

(kglm3)  (%) 

5E 

(Si) 

si 

(Si) 

a 

(m3 Is) 

n 

k 

(m2) 

I'. 

(m3) 

1 

84,000 

2.15X1 0'3 

470 

51.9 

5.7 

3.0 

1.1X1  O'7 

439 

2.47X1  O'8 

3.2X1  O'3 

2 

99,300 

6326x1 04 

490 

483 

3.8 

3.0 

1.6X1  O'8 

3.07 

2.55X1  O'9 

9.41  X104 

3 

162,300 

5.049X1 0'3 

470 

52.8 

7.1 

5.0 

2X1 0'7 

3.29 

4.96X10'™ 

6.19X10'3 

4 

327,900 

3.28X1 0'3 

430 

573 

6.2 

5.0 

53X10-8 

2.85 

2.1 3X1  O'70 

3.92X1  O'3 

5 

326320 

3.49X1  O'3 

430 

57.5 

6.2 

5.0 

4X1  O'8 

233 

6.58X10'11 

4.22X  10'3 

6 

349,800 

4.86X1 0'3 

510 

48.4 

8.1 

6.0 

2.25  XI 0'7 

3.79 

1.41X10'9 

5.92X1  O'3 

7 

255,300 

4.09X1 0'3 

570 

48.7  8.2 

6.0 

3X1  O'7 

4.31 

9. 75X1  O'9 

5. 19X1  O'3 

8 

514,800 

4.42X1  O'3 

499 

48.5 

5.7 

5.5 

1.23X1  O'6 

2.38 

6.57X10'10 

4.49X1  O'3 

9 

265,000 

1.87X10'3 

520 

46.2 

5.3 

5.0 

3.05  XI  O'7 

2.43 

8.72X1  O'10 

2X1  O'3 

10 

249,840 

3.1  XI 0'3 

490 

51.3 

8.1 

8.0 

6.2X1  O'7 

2.16 

4.58X10'™ 

3.15X10'3 

11 

221,400 

3.31  XI  O'3 

500 

50.4 

8.06 

5.0 

4.18X1  O'7 

3.39 

2.5  74X1  O'9 

4.89X1  O'3 

i 


limits  of  the  unsaturated  flow  theory  used  here  should 
be  examined  by  labeling  the  water  with  tracers  and 
doing  experiments  under  conditions  where  the  snow 
properties  can  be  observed  to  change  throughout  the 
duration  of  the  experiments. 
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Figure  A 7.  The  measured  water  discharge  points  are  shown  against  time  along  with  cal- 
culated profiles  of  discharge  versus  time  for  various  values  of  n for  Experiment  no.  7. 


Figure  A 8.  The  measured  water  discharge  points  are  shown  against  time  along  with  cal- 
culated profiles  of  discharge  versus  time  for  various  values  of  n for  Experiment  no.  9. 
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